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We study the formation of Majorana states in quasi-ID and 2D square lattices with open boundary 
conditions, with general anisotropic Rashba coupling, in the presence of an applied Zeeman held and 
in the proximity of a superconductor. For systems in which the length of the system is very large 
(quasi-ID) we calculate analytically the exact topological invariant, and we hnd a rich corresponding 
phase diagram which is strongly dependent on the width of the system. We compare our results 
with previous results based on a few-band approximation. We also investigate numerically open 2D 
systems of hnite length in both directions. We use the recently introduced generalized Majorana 
polarization, which can locally evaluate the Majorana character of a given state. We hnd that the 
formation of Majoranas depends strongly on the geometry of the system: for a very elongated wire 
the hnite-size numerical phase diagram reproduces the analytical phase diagram for inhnite systems, 
while if the length and the width are comparable no Majorana states can form, however, one can 
show the formation of “quasi-Majorana” states that have a local Majorana character, but no global 
Majorana symmetry. 
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I. INTRODUCTION 

Many theoretical works have predicted the existence 
of Majorana fermions in various solid state systems,^ 
for example in one-dimensional semiconducting wires 
with either strong spin-orbit coupling in the presence 
of an uniform Zeeman magnetic held,^“^ or without 
spin orbit but with an inhomogenous Zeeman field, 
as well as in various two-dimensional systems and 
arrays.Though there have already been several 
promising experiments^^”^^ their interpretation is con¬ 
troversial and no clear-cut detection of Majoranas has 
so far been attained. An important aspect to under¬ 
stand in order to connect with the experimental mea¬ 
surements is the effect of the geometry of a system in 
the formation of Majorana states. Most previous works 
dealing with the formation of Majoranas have focused 
on either one-dimensional finite-size systems or on two- 
dimensional nanoribbons with periodic boundary condi¬ 
tions, with a few exceptions,some of them focusing on 
multi-channel one-dimensional wires. 

Here we are interested in the effect of the geometry of 
a system on the conditions to form Majorana states in 
quasi-ID and 2D systems. We focus on two main ques¬ 
tions: what are the conditions for the formation of Majo¬ 
ranas 1) in quasi-ID infinite systems with both longitu¬ 
dinal and transverse Rashba coupling, and 2) in systems 
with open boundary conditions in both directions, also in 
the presence of both longitudinal and transverse Rashba. 

To answer the first question, previous works have fo¬ 
cused on using multi-band approximations,or 
on perturbing a simpler underlying model,both ap¬ 
proaches limiting though the applicability of the results 
to a small region of parameter space. Here we use a differ¬ 


ent technique which allows us to go beyond these approx¬ 
imation and calculate exactly the topological invariant 
for a system with arbitrary longitudinal and transversal 
hopping and spin-orbit coupling, using a fully analyti¬ 
cal technique similar to the one presented in Ref. 43. 
This calculation allows us to determine the topological 
phase diagram for such systems which is very rich and 
depends strongly on the width of the wire. Moreover, 
we note that, in contrast to the single channel model the 
strength of the spin-orbit coupling now enters explicitly 
into the form of the invariant and affects the form of the 
phase diagram. We compare our exact results with re¬ 
sults obtained using an approximation of our system to 
a multi-band one-dimensional system, which is however 
valid only for a small region in the parameter space. We 
can show that in this region of the parameter space the 
exact topological phase diagram is recovered accurately 
by the multi-band approximation, however our theory 
is much more general and can describe the topological 
phase diagram for the entire parameter space. 

The second question that we address, the condition 
for the formation of Majoranas in fully open systems of 
arbitrary size, has, to our knowledge, not been previ¬ 
ously addressed quantitatively for realistic models. Pre¬ 
vious studies have focused on purely one-dimensional 
systems, two-dimensional systems with periodic bound¬ 
ary conditions (i.e. infinite nanoribbons),as well as 
open systems in which the aspect ratios are very large 
(quasi-ID) by using certain approximations, as men¬ 
tioned above.However, the case of the fully 
open systems of arbitrary size has not been previously 
addressed, and the analysis of such system is much more 
complex, as we expect that in the presence of fully open 
boundary conditions the Majorana states would connect 
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and hybridize, and thus get destroyed. This phenomenon 
has already been seen in the quantum Hall effect.More¬ 
over, such analysis is more relevant than the previous 
studies for the actual experimental situations'^ , since 
the wires studied in the transport experiments are neither 
purely one-dimensional, nor infinite ribbons, but open 
quasi-2D or quasi-3D systems with a finite transverse 
width for which transversal hopping, as well as transver¬ 
sal Rashba coupling are important. 

To explore this complex problem we vary the ratio 
between the length and the width of our system and 
we analyze the low-energy states numerically, using the 
MathQ code^^. A crucial tool in our analysis is the Ma- 
jorana polarization (MP) whose general definition was 
recently introduced in Ref. 48; this plays the role of a 
local measure of the Major ana character of a given state. 
Previous studies have focused only on the local density 
of states, which carries no direct information about the 
Majorana nature of a state.Using the MP we find that 
the formation of Majorana states in fully open systems is 
strongly affected by the geometry of the system. Thus, 
for asymmetric elongated systems, i.e. the quasi-lD limit, 
we recover numerically the phase diagrams that we had 
obtained also analytically by calculating the topologi¬ 
cal invariant. When the system is deformed towards a 
quasi-square system, the Majorana states are gradually 
destroyed for larger and larger regions in the parame¬ 
ter space, and for a square system Majorana states can¬ 
not form, however low-energy quasi-Majorana states ex¬ 
ist. Such states have locally an almost perfect Majorana 
character, however the Majorana polarization vector is 
not aligned but varies spatially and thus a full Majorana 
state cannot form. Such states have also been observed 
in the spinless Kitaev model.The parameter range in 
which such quasi-Majorana states form corresponds sur¬ 
prisingly enough to the 2D bulk topological phase pre¬ 
dicted for this system.Note that these states are 
similar to the ‘chiral’ Majorana states,which are 
also Majorana-like low-(but non-zero)-energy states aris¬ 
ing on the internal and external boundaries of a finite 
size system with a hole. These ‘chiral’ Majorana states 
would become exact Major anas when a magnetic fiux is 
threaded through the system. The relationship between 
the quasi-Majorana states and the chiral Majorana states 
is the topic of a separate investigation. 

We should note that the two approaches described 
above for the two questions addressed in this work are 
very general and can be generalized to predict the con¬ 
ditions for Majorana formation for systems of arbitrary 
size, dimension and lattice structure, as well as for any 
type of hopping and spin-orbit couplings. While we focus 
here on a rather simple tight-binding model applying to a 
two-dimensional square lattice, we should stress that our 
approach is fully generalizable to more complex tight- 
binding models that can map more accurately realistic 
systems such as In As and InSb wires; a full numerical 
analysis of such more complex models is envisageable in 
the future. 


This paper is organized as follows. In Sec. II we in¬ 
troduce the general model. In Sec. Ill we calculate an¬ 
alytically the topological invariant for a generic infinite 
quasi-ID systems. In Sec. IV we present the numerical 
phase diagrams obtained for finite-size systems of differ¬ 
ent geometries using the Majorana polarization, and we 
show that in the limit of very elongated wires this coin¬ 
cides with the phase diagram computed analytically in 
Sec. III. In Sec. V we look at the signatures of the phase 
transitions in experimentally measurable quantities, such 
as the density of states. In Sec. VI we explore the connec¬ 
tion with the previously studied multi-band models and 
we show how we can recover the same results for specific 
regions in the parameter space. We conclude in Sec. VII. 
More detailed examples are presented in the Appendices, 
such as the effects of the asymmetry between the longi¬ 
tudinal and transverse Rashba, as well as the transition 
between open and periodic boundary conditions. 


II. MODEL 

We first introduce the Hamiltonian which under¬ 
lies our analysis. We use the Nambu basis: Tt = 

{ct^, —c^q}, where annihilates (creates) a 

particle of spin a at site r = {i,j) in a square lattice. The 
corresponding wavefunction is 

We also use r to denote the Pauli matrices in the particle- 
hole subspace and a for the Pauli matrices in the spin 
subspace. We consider the two dimensional model Hamil¬ 
tonian 






+4']^ {-tx - ictxcr^) + H.c. 

+4't {-ty + + H.c. 


( 1 ) 


tx,y is the nearest neighbour hopping, /i the chemical 
potential, B the magnetic field, A the proximity in¬ 
duced s-wave superconducting pairing, and ax^y are the 
strength of the longitudinal and transversal Rashba spin- 
orbit couplings. We set = t = 1 and h = 1 through¬ 
out. Unless explicitly stated we also set ty = tx = 1. 
Here x and y are the unit vectors for the x and y direc¬ 
tions respectively, and the lattice spacing is a = 1. We 
will focus exclusively on square lattices. We also define 
the particle hole operator C = e^^ where K is 

the complex-conjugation operator, and is an arbitrary 
phase. The Hamiltonian anti-commutes with this opera¬ 
tor, {C,H} = 0, and = 1. 

The size of the system to be analyzed is Nx x Ny. For 
the quasi-lD systems Nx ^ Ny. The magnetic field is 
taken to be along z. For the 2D open systems we consider 
here a magnetic field perpendicular to the plane and we 
neglect the orbital effects; the orientation of the magnetic 
field may affect the results but we do not focus on these 
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FIG. 1. (Color online) Bulk topological phase diagram for (a) 
a strictly ID - single channel, wire and (b) a square lattice, 
where A = 0.4t. White is the topologically trivial phase, light 
red is the bulk topologically non-trivial phase where a pair of 
Majorana edge states form at the boundaries. The black lines 
show where the bulk gap closes. 


effects in this manuscript, but it is the focus of a separate 
investigation. 

Two limits are already well understood. If = 1 we 
are back to a strictly ID wire. If both the x and y direc¬ 
tions are bulk, one has the standard 2D case.^’^’^^’^^’^^ 
The topological phase diagram of both of these cases is 
already known, see Fig. 1. The 2D array is in the topo¬ 
logical class D described by a Chern number, a Z invari¬ 
ant. The ID case is in a higher symmetry class BDI^^ 
due to the absence of ay and is also described by a Z 
invariant, a winding number. The introduction of ay in 
the quasi-ID breaks a so-called time-reversal symmetry 
putting it into the D class which in ID has a Z 2 invariant, 
calculated from the Pfafiian at high symmetry points. 

As we have only nearest neighbor terms, the invariant in 
ID BDI, i.e. the winding number, only takes the values 
u = 0, ±1. One can also define the parity of a Z invari¬ 
ant, 5 = e^^^ = ±1. This is generally sufficient for our 
purposes as the quasi ID systems in the D class have a 
Z 2 invariant and in this case S and u have the same infor¬ 
mation. For the Z invariant limited to u = 0, ±1, the lost 
information is small and not crucial to our arguments. 

In 2D the Chern number can take several values, 
though we limit ourselves to plotting only the parity of 
the invariant in Fig. 1. For both the ID and 2D sys¬ 
tems the topological phase diagram is independent of the 
strength of the Rashba coupling, except for the condition 
that it must be non-zero. In practice, very small Rashba 
couplings would require very large systems to see clearly 
defined MBS. As we shall see in the following section, 
the quasi-ID systems have a much richer phase diagram, 
which also depends explicitly on the strength of ay. An¬ 
other example where such dependence on the spin-orbit 
strength is a hexagonal lattice with a similar Hamiltonian 
to the one described in Eq. ( 1 ). 28 , 43,53 

Previous models were limited to either an approximate 
low band theory^'^’^^’^^^'^^ or neglecting the role of spin- 
orbit coupling all together.We go beyond these approx¬ 
imations here by calculating the topological invariant of 


the Hamiltonian given in Eq. (1) exactly. This allows a 
much more accurate depiction of the phase diagram, in 
particular for large magnetic fields and large spin-orbit 
coupling where the approximate models do not apply. 
The experiments on nanowires deposited on the top of 
superconducting substrates^^ are more accurately cap¬ 
tured by this model, as it is not in general possible to 
neglect the role of transverse spin-orbit coupling. 

III. TOPOLOGICAL PHASE DIAGRAMS FOR 
QUASI-ID WIRES 

We will start by computing analytically the exact 
phase diagram for quasi ID systems (finite width and in¬ 
finite length, with open boundary conditions in the trans¬ 
verse direction). For a single wire Ny = the topologi¬ 
cal phase diagram is presented in Fig. 1. When Ny > 1 
the topological phase diagram can be obtained from a 
calculation of the topological invariant of the system, as 
described in what follows. 

It is often possible to find a Z 2 invariant for a topologi¬ 
cal superconductor with particle-hole symmetry by map¬ 
ping the problem to the parity at the time reversal in¬ 
variant momenta (Ti).^^ This is equivalent to calculating 
the Pfafiian of the Hamiltonian^’^^ but proves more con¬ 
venient for the purposes here. Parity in this case refers 
to an operation which commutes with the Hamiltonian 
at these points, [P,H{Ti)] = 0, and anti-commutes with 
the particle-hole operator, {P, C} = 0. The eigenstates of 
the Hamiltonian have parity eigenvalues ±1 at the TRI 
momenta. It follows that there is a representation in 
which the Hamiltonian is block diagonal at F^ with two 
blocks, one with parity 1 and the other with parity —1. 
This is the most convenient basis for calculating the in¬ 
variant, and for a ID wire or 2D square lattice a trivial 
re-ordering of the Nambu basis suffices. In the system 
that we will consider here the problem can be reduced to 
finding this appropriate basis so that the Hamiltonian is 
block diagonal at F^ . 

Firstly we note that if ay = 0 then the calculation 
becomes straightforward and standard, see App. B or 
Refs. 22-25, and 28. This is because the problem can 
be decomposed into the transverse momentum channels 
labelled by the quantum number ky. For each ky the 
Hamiltonian at the time reversal invariant (TRI) mo¬ 
menta, Fi = 0 and F 2 = tt, is trivially block diagonaliz- 
able, exactly as for a ID wire, in such a way as to define 
a parity operator which both commutes with the Hamil¬ 
tonian at the TRI momenta and anti-commutes with the 
particle-hole operator C. All that is required is a re¬ 
ordering of the Nambu basis to {ct^, ct^, —c^^}. In 
this basis the parity operator is simply 

=(0 -i) ^ p) 

Here is the n x n identity matrix. This is sufficient for 
calculating the topological invariant. 
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For the quasi ID system with transversal Rashba con¬ 
sidered here the Hamiltonian is no longer block diagonal- 
izable in any trivial way at the TRI momenta. Nonethe¬ 
less it is possible to find a suitable parity operator, or 
equivalently a suitable transformation on the Hamilto¬ 
nian, by using the transverse spatial, spin and particle- 
hole subspaces. 

After a Fourier transform along the wire the Hamilto¬ 
nian can be written as with 


n{k) = 


/ fk ^ B 
C) 
-A 

V 0 


£^k 0 \ 

fk-B 0 -A 

0 B - /I, 

-A C*_, -fU-B) 


( 3 ) 


where the entries are themselves matrices for the trans¬ 
verse spatial direction y. The terms B and A are diagonal 
in this space and 
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= u la L.k -\a ... ( 5 ) 


Finally /(A:) = —‘Itx cos[fc] — and £k = —2iax sin[A:]. 

The transformation that allows for the Hamiltonian to 
be written in block diagonal form at the TRI momenta 
is 'H(k) = W'H{k)U where 


U = 




+ 


cr"T" 


- , 


(6) 


where we introduced an Ny x Ny matrix whose elements 
are given by = 5n,Ny+i-n'- Then 


n^i) 


mti) 0 \ 

V 0 • 


In this basis the parity is 


PNy = 



0 ^ 

-l2ivd ’ 


( 7 ) 


(8) 


and we have both [P/v^,7i(f^)] = 0 and {P^y^C} = 0. 
We note that also [CM] = 0- 

The topological invariant is then calculated in the 
usual way^ 


( 9 ) 
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FIG. 2. (Color online) Topological phase diagram as a func¬ 
tion of B and y for a ladders with (a) Ny = 2, (b) Ny = 3, 
(c) Ny = 15, and (d) Ny = 50 with A = 0.4. Light red is 
the topologically non-trivial phase and white is topologically 
trivial. For (a,c) Occ = = 0.2 and for (b,d) 

For (a,b) gap closing lines belong to the two different TRI 
momenta are shown as dashed and full lines . In (c) and (d) 
the gap closing lines for the 2D bulk system are shown for 
comparison. 


When ^ = — 1 there is band inversion, i.e. the parity 
switches between the TRI momenta an odd number of 
times and the system is topologically non-trivial. For 
6 = 1 the system is topologically trivial, see Fig. 2 for 
some examples. 

The full analytical expressions become quickly cumber¬ 
some as Ny is increased. Full expressions for Ny = 2,3 
are give in App. A. One key point to note is that be¬ 
cause of the transverse spin orbit coupling the expressions 
detP(ri^ 2 ) do not factorize further into simple expres¬ 
sions for each wire. 

The two wire system (two-leg ladder) presented in 
Fig. 2(a) has a peculiar property: close to the point 
B = 2t and /i = 0 inside the topologically trivial phase, 
two pairs of MBS can form. These states are not how¬ 
ever topologically protected and can be destroyed by in¬ 
troducing disorder.The zero overlap of these MBS is 
sometimes referred to as a “hidden symmetry”.^^ 

In contrast to a wire or square lattice, but similar to 
a hexagonal lattice,the topological phase diagram 
depends explicitly on all of the parameters including the 
Rashba spin-orbit coupling. Thus we note that the num¬ 
ber of the topological regions in the area close to /i = 2, 
and B <2 (between the ID and 2D transition lines) de¬ 
pends on both the width of the system, as well as on the 
value of the Rashba coupling, and it strongly reduced 


6 = sgndet H(f i) det H(f 2 ). 
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for large Rashba values as well as for a large Ny. One 
must thus note that switching on the spin-orbit coupling 
while keeping the other parameters constant, can also 
drive the system through a topological phase transition. 
Also, while in this section we focus mostly on the case 
of isotropic coupling, anisotropic systems in which the 
spin-orbit interaction and hopping strengths may not be 
symmetric between the x and y directions are explored 
in App. C. We note that for finite when ^ 0 
we recover the phase diagram for an array of wires, see 
Fig. 8. In this limit the nanowire can host many topologi¬ 
cally protected Majorana bound states. Switching on 
which breaks a time reversal symmetry of the problem, 
gaps out all but at most two of the MBS. Naturally in the 
weak coupling limit it is possible to infer the topological 
phase diagram from that for ay = 0.^^ 

Note the evolution from the phase diagram for a sin¬ 
gle wire shown in Fig. 1 to the complex pattern that we 
obtain for a large Ny. It is important to stress here the 
difference between the phase diagram for a system with 
open boundary condition in which Ny is very large (such 
as Fig. 2(d) ), and the bulk phase diagram for the 2D 
system depicted in Fig. 1(b). The latter corresponds to 
the topological phase diagram describing the formation 
of Majorana edge state for a system with infinite Ny and 
in which periodic boundary conditions are imposed in the 
y direction. The main difference consists in the fact that 
in the system with open boundary conditions the phase 
diagram remains discreet, and the trivial and topological 
phases alternate with y and B with periodicities which 
become smaller and smaller with increasing the width of 
the system. The source of this discreteness is the fact 
that, no matter the size of the system, it is still an open 
system with a single boundary. The open system, invari¬ 
ant of dimension, is fundamentally different from the in¬ 
finite nanoribbon (periodic system), and we believe that 
this difference has not been thoroughly explored, espe¬ 
cially in connection with experimental systems such as 
quasi ID wires, which should rather be modeled as finite- 
width open systems. A more detailed discussion of the 
difference between fully open systems (flakes) and sys¬ 
tems with periodic boundary conditions (tubes) as well 
as the transition between the two cases, is presented in 
Appendix C. 

Note also that the existence of the discrete steps in 
the phase diagram (known also as the ‘even-odd’ effect) 
has been also predicted for the multi-band wires:the 
relation between our exact calculation valid for wires of 
arbitrary widths and in all the parameter space and the 
multi-band approximation is explored in detail in Sec. VI. 


IV. 2D OPEN SYSTEM 

In what follows we consider fully open systems of finite 
dimensions in both directions, and we analyze their prop¬ 
erties numerically by looking at the lowest energies eigen¬ 
states obtained using the MathQ code.^^ If Ny ^ Nx 


one expects to recover the phase diagrams obtained an¬ 
alytically in the previous section. The interesting ques¬ 
tion is what happens when the transversal and longitu¬ 
dinal directions become comparable, in which situation 
we should expect that the system has effectively a single 
boundary and the Majorana states on the edges would 
hybridize and destroy each other. The manner in which 
this happens, as well as the persistence of some topologi¬ 
cal properties in the low-energy edge states of a finite-size 
open system are investigated using the generalized Ma¬ 
jorana polarization introduced in Ref. 48. 


A. Phase diagram deduced from the Majorana 
polarization 

While for an infinite system the condition to be in a 
topological state can be determined using the topological 
invariant, such criterion cannot be applied to finite-size 
systems. This can happen for example because finite-size 
effects can destroy the Majorana states even in regions 
in the phase space that are characterized as topological 
from a bulk perspective. Thus one needs to apply a dif¬ 
ferent criterion for such systems to test the formation of 
Majoranas. Moreover, solely the existence of zero-energy 
states is not enough since many non-Majorana states, 
such as impurity states, can form at zero energy; also 
many infinitesimally-small-energy states can arise which 
are hard to distinguish from purely zero-energy states. 
To test if such states are Majorana or not we are apply¬ 
ing here the criterion of the generalized Majorana polar¬ 
ization (MP)."^^ 

The MP is a local measure characterizing a specific 
wavefunction which can tell if such state is a Majorana 
or not. While the topological invariant is a character¬ 
istic of the ‘Hamiltonian’ of a given system, the MP is 
a characteristic of the ‘eigenstates’ of the Hamiltonian, 
as calculated for a specific configuration, for a system 
of a given size and with particular boundary conditions. 
Thus, in the same way that the local density of states 
(LDOS) can provide space-dependent information about 
the formation of edge states for a given configuration of 
the system, information which is not solely dependent on 
the Hamiltonian, and cannot be extracted directly from 
it, the MP can tell about the formation of Majoranas in 
a given system by analyzing its eigenvalues and eigen¬ 
functions, and their dependence on position. The MP 
is thus an equivalent of the LDOS measuring the Majo¬ 
rana density at a given position in space. However, while 
the LDOS is a scalar quantity, the MP is a vector in 
the complex plane (we can see it as a Majorana pseudo¬ 
spin), whose magnitude describes the Majorana density, 
i.e. the electron-hole overlap, and its phase the relative 
phase between the electron and hole components. 

The way in which the MP criterion works is as 
follows:"^^ a Majorana state is an eigenstate of the particle 
hole operator. Therefore a Majorana-like state localized 
inside a spatial region IZ must satisfy (7 = 1 where C is 
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the magnitude of the integral of the Major ana polariza¬ 
tion vector over the spatial region R. 




( 10 ) 


f is the projection onto site r, and the local MP is simply 
the expectation value of the local particle-hole transfor¬ 
mation: 


(^1 C f 1^) = -2 UUf^crVrcr ■ (H) 

CT 

We split our system in half such that IZ is half of the sys¬ 
tem (divided usually along the longer length). The Ma- 
jorana polarization also allows us to distinguish between 
zero energy states which are Majorana bound states and 
those which are not as C = 1 is both a necessary and 
sufficient condition for a state being a Majorana state. 
By numerically solving the tight-binding Hamiltonian, 
Eq. (1), and plotting C as a function of the parameters we 
can accurately recover the appropriate topological phase 
diagram. The topological phases are characterized by 
(7 = 1, while the non-topological ones by (7 = 0. This 
is demonstrated in Fig. 3 for a few systems of different 
sizes. 

The MP thus allows an accurate characterization of the 
topological phase even in cases where a direct calculation 
of the topological index is not feasible. Moreover it does 
not require very large system sizes to be highly accurate, 
as for example a map of the lowest energy state would. In 
addition, as we will shortly demonstrate, the MP gives 
local information about the formation and behavior of 
the Majorana bound states which can not be gained form 
a consideration of the system’s bulk properties. Thus we 
should note that, while if a Majorana state exists the 
system is surely topological, the reverse is not necessary 
true, as a bulk topological system may have no Majorana 
solutions in certain finite-size configurations. 

For a system in which the length is much larger than 
the width Ny Nx such as the one depicted in Fig. 3(a) 
we would expect to recover the phase diagram calculated 
analytically in the previous section, and this is indeed the 
case, the topological regions predicted by the MP crite¬ 
rion (red in Fig. 3(a)) correspond exactly to the topolog¬ 
ical regions (light red) in Fig. 2(b). For a wider system 
such as the one depicted in Fig. 3(c) we can see that the 
phase diagram calculated analytically (Fig. 2(d)) is re¬ 
covered very well at large values of B, however, at small 
values of B, the value of C is neither 0, nor 1, but it shows 
an intermediate value (depicted in yellow). For fully 
square systems (Fig. 3(b)) the entire topological phase 
is characterized now by (7 ^ 0.7. We should note that all 
the states corresponding to the regions with (7 < 1 have 
non-zero but very small energies, thus making it harder 
to interpret them as Major anas or not based solely on 
the energy criterion. 

To understand these states we focus on the local struc¬ 
ture of the MP defined in Eq. (11) as a vector in the com¬ 
plex plane. A Majorana state must have an integral of 


(a) (b) 
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FIG. 3. (Color online) The phase diagram as a function of 
B and for open systems of (a) 3 x 201, (b) 101 x 101, and 
(c) 51 X 201. This has been calculated using the PH expec¬ 
tation value of the lowest energy state summed over half of 
the system C (see main text). In all the examples A = 0.4t 
and a = 0.5t. The analytical topological phase boundaries for 
quasi ID (solid lines) and 2D (dashed lines) are also plotted. 


the MP of 1 over a spatial region 7^, thus it must exhibit 
an uniform phase inside this region, equivalently the MP 
local vectors need to be aligned inside IZ (‘ferromagnetic’ 
MP structure). In Figs. 4-5 we plot the MP for a variety 
of different low energy states. In Fig. 4 we plot the MP 
vector for a 51 x 201 system with A = 0.3t and a = 0.5t. 
Fig. 4(a) corresponds to the red (topological) phase in 
Fig. 3(c) ((7 = 1), with /x = 3.5t and B = 2.2t. We should 
expect to have two MBS confined at the two narrow ends 
of the ribbon, and we note that this is indeed the case, 
we can see the formation of two ‘ferromagnetic’ states 
localized at the two ends of the wire, with opposite MP. 
Fig. 4(b) corresponds to a non-topological (blue phase 
in Fig. 3(c)) phase, with (7 = 0. We take /i = 3.5t and 
B = 2.3t, and the corresponding state, while being lo¬ 
calized on the edges, is indeed non-Majorana, which can 
be seen from the local fluctuations of the direction of the 
MP vector from site to site, which makes its integral over 
half of the wire zero. In Fig. 4(c) we focus on the more 
puzzling (7 < 1 phase, denoted in yellow in Fig. 3(c). We 
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FIG. 4. (Color online) MP for the lowest energy states for 
open system of 51 x 201 with A = 0.3t and a = 0.5t. We plot 
(a) a MBS fi — 3.5t and B = 2.2t (corresponding to a red 
phase in Fig. 3(c)), (b) a ID edge state = 3.5t and B = 2.3t 
(blue phase in Fig. 3(c)), (c) a 2D edge state = 3.5t and 
B — t (yellow phase in Fig. 3(c)), and (d) a bulk state /x = 0.5t 
and B = 5t (blue phase in Fig. 3(c)). 


take /i = 3.5t and B = t. We note that the corresponding 
lowest energy state is also an edge state, extending over 
the entire contour of the system. Moreover, while the 
state is not ‘ferromagnetic’, a small uniform variation is 
observed from site to site, making the total integral over 
R finite, but not equal to 1. Finally in Fig. 4(d) we show 
a state with (7 = 0 for /x = 0.5t and B = bt, which corre¬ 
sponds to the blue region outside the region delimited by 
the bulk 2D topological lines, this actually correspond to 
a bulk state for which the value of the MP is small even 
locally. 

It appears that the existence of the ‘yellow’ phase is 
due to the shape of the system, i.e to the fact that the 
length and the width become comparable. If the length is 
increased the yellow phase is diminishing and the phase 
diagram converges towards the phase diagram calculated 
analytically in the previous section. This observation has 
been made also for the Kitaev model.However, we want 
to stress that in real experimental systems the existence 
of the intermediate phase with (7 < 1 is possible, and it 
would correspond to the formation of low-energy subgap 
states which have locally a full Major ana character, but 
that are non-Majoranas, due to their global lack of sym¬ 
metry. It would be interesting to explore the connection 


101 


51 


1 



1 51 101 


X 



FIG. 5. (Golor online) MP for lowest energy state for open 
system of 101 x 101 with A = 0.3t and a = 0.5t. (a) fi = 4t 
and B — and (b) /x = 2t and B = 3t, both in the 2D 
topologically non-trivial phase. 


between these states and the ’chiral’ Majorana states de¬ 
scribed in Refs. 19 and 50. Another interesting question 
would be if such quasi-Majorana states have any topo¬ 
logical characteristics, such as non-Abelian statistics, or 
atypical braiding properties, and if such states can be 
useful for example for quantum computation in a simi¬ 
lar manner as the Majorana states. Note that a state 
with a (7 = 1 is a Majorana, and certainly has non- 
Abelian statistics. It is not at all clear that the statistics 
of the quasi-Majorana states would be Abelian. These 
are peculiar states in-between fermionic and Majorana- 
like, and there exist to this point no work to investigate 
their statistics. Moreover there seems to be a fundamen¬ 
tal difference between the properly fermionic states which 
have a uniform local zero MP. The quasi-Major anas have 
their weight divided into a few different regions, roughly 
disconnected from each other, each with C different from 
zero. Manipulating such a state in one part of the system 
is not the same as manipulating the fermionic state, and 
the effect of such manipulation is also a very interesting 
question. Such questions are however well beyond the 
scope of the present work. 


An extreme situation is depicted in Fig. 5(a,b), in 
which we plot the MP for two quasi-Majorana states for 
a square system. Fig. 5(a) depicts the MP for a set of 
parameters inside the 2D bulk topological phase (as de¬ 
scribed in Fig. 1(b)), while Fig. 5(b) describes a set of pa¬ 
rameters inside the ID topological phase (as described in 
Fig. 1(a)), that would become a Majorana end state if one 
direction of the system was sufficiently increased. Inter¬ 
estingly enough these two states have different MP char¬ 
acteristics, the ‘ID’ one being localized roughly at the 
corners, while the ‘2D’ one along the contour of the sys¬ 
tem, allowing one to distinguish between the two. Note 
that for a fully square system the corner MP vectors are 
oriented at roughly 90 degrees angle with respect with 
each other, giving rise to an integral of (7 = V2/2 « 0.7, 
which corresponds well to the ‘yellow’ in Fig. 3(b). 

























V. LOCAL DENSITY OF STATES 


It would be interesting to see the effects of the dis¬ 
crete topological transitions in the phase diagram of the 
open systems on directly observable quantities such as 
the local density of states (LDOS) measurable in STM 
or tunneling experiments. This has been explored for 
systems with a small number of bands in certain regions 
of the parameter space,but in order to make 
contact with the experiments we need to explore systems 
of different sizes and geometries. In what follows we fo¬ 
cus on both quasi-ID systems, as well as on systems in 
which the width becomes comparable to the length of the 
wire. The LDOS at one end of the wire is defined as 

/3(e) = + T^)fedge| V’n) ^ 2 

r is a broadening energy scale, typically taken to be 
slightly larger than the mean level spacing. To avoid 
effects due to local fluctuations we define the edge LDOS 
as an average over all sites on the first 5 atomic rows 
away from the edge. 

In Fig. 6 the LDOS at one end of the quasi-lD wire is 
shown for two examples, (a) 15 x 201 sites and (b) 51 x 201 
sites. In both examples A = 0.4t, /i = 4t, and a = 
0.5t. The bulk gap closing is only faintly visible visible 
in Fig. 6(a) due to the small weight of the bulk states at 
the edges of the system. We note the formation of MBS 
edge states inside the parameter regions predicted by the 
topological phase diagram. In addition, we note that 
even in the topologically trivial phases, the lowest energy 
states are also edge states, even if split in energy and 
non-Majoranas (see also the characterization of such a 
‘blue’ state in the previous section). Note that these edge 
states can be understood as MBS of a BDI model split 
by the transverse Rashba, see App. B. As the strength of 
the Rashba coupling is relatively small in the nanowires 
currently used^^ this splitting could be experimentally 
very small. 


VI. EFFECTIVE MULTI-BAND THEORY 



0 1 2 3 0 


p(e) (a 1) 



FIG. 6. (Color online) The local density of states at the edge 
of quasi-ID wires as a function of energy and B with A = 0.4t, 
/i = 4t, and a = 0.5t. Panel (a) is for Ny = 15 and (b) is for 
Ny = 51 , in all cases Nx = 201. 


~i^xx + with m being the effective mass. Dis¬ 

cretizing on a lattice scale a we find the lattice Hamilto¬ 
nian Eq. (1) with a corresponding hopping t = l/(2ma^) 
and a scaled chemical potential of dji = /i — /io, and 
/io = 4t. Thus J/i = 0 corresponds to the point where we 
can define our effective model which is the bottom of the 
band. 

To begin we substitute 






2 . nny 

sm ■ 


n=l 


Ny “h 1 Ny + 1 


(13) 


into Eq. (1). Then we have 




QID 


= E 


(fnT^ - - At® + Bcr®) ?/)„ 


Eor certain ranges of parameters the phase diagram 
for a fully open system can be calculated using a multi¬ 
band approximation. Such an approximation is valid if 
the width Ny is large, respecting however the quasi ID 
condition, N^ ^ 19 , 34 , 38,40 However it is clear that 

much of the rich physics that we can find exactly for 
system of arbitrary width and length cannot be recovered 
using this approximation. In what follows we show how 
one can recover the appropriate multi-band limit for our 
model, based on a truncation of the active bands in the 
model, and we show how it compares to the exact result 
in terms of the regions of validity. 

Let us start by relating our tight-binding model to 
a continuum model with a quadratic kinetic energy: 


+'>Pix {-t - T^i)n,x+1 + H.C. 

“ E '^ix^Nxr'"r^Mnn''^n'x , 

x^ 

nn' 

where = —2tcos[kn], = 7rn/{Ny + 1) and 


(14) 


Afnn' — 


0 

sin[(/cr,+/c^/)/2] 

sin[(/c^-/c^,)/2] 


if n + n' is even. 


sin[(kn — k^/ )/ 2 ] 
sin[(/cri+/c^/)/2] 


otherwise. 


(15) 

By truncating the transverse modes at some particular n* 
we obtain an approximate theory consisting of just 4n* 
bands rather than 4V^. We will refer to n* as the num¬ 
ber of channels, and each channel carries both spin and 
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particle-hole degrees of freedom. The energy difference 
= Cn* — Cl gives an energy scale for the multi-band 
theory. We need to impose B^A^a £ri*^ where for 
large wire widths {Ny ^ 1) we can write 


" L 


{n*f - 1 


(16) 


Naturally the mixing between the bands with n < n* 
and n > rf must be negligible compared to , i.e. the 
Rashba coupling must be small. 

In order to calculate the topological invariant we pro¬ 
ceed exactly in the same way as for the exact calcula¬ 
tions presented in Sec. II, and we start by performing a 
Fourier transform along the bulk direction of the wire. 
The Hamiltonian becomes H = 
with Unn'ik) = 'Hn(k)5nn' + ^nd 



/ fkn + B 


-A 

0 

Kik) = 

r) 

-A 

fkn 

0 

e 

aT 

o 

1 

-A 

^-k 


\ 0 

-A 

^ * 

^—kn 

J —k,n 

where 

^ = —2tcos[k] - 

- 2tcos[kn] — 

6fi and 


—2\ax sin[/c]. The band mixing terms are 


\ 


-bJ 

(17) 

= 


qiy _ 


/ 0 Mnn' 0 0 \ 

2iay [ Mnn' 0 0 0 I 

Ny + l 0 0 0 Mnn' 

Vo 0 Mnn' 0 / 


(18) 


The transformation that allows for the Hamiltonian to 
be written in block diagonal form at the TRI momenta 


is Brin' ik) = UlUnn' {k)Un' where 


1 — (J^T^ 

Un = - ^ + 


CT^T" — 


(19) 


Following the same procedure as the one presented in 
Sec. II recovers the phase diagram more accurately. In 
Fig. 7(c) we compare the approximate phase diagrams 
for systems with different numbers of wires, Ny = 100 
and Ny = 101, obtained using the same number of chan¬ 
nels, n* = 10. The sole difference between the two phase 
diagrams are the small extra topological regions denoted 
in pink for Ny = 101 which are barely visible, confirming 
that the multi-band theory is not very sensitive to Ny in 
its region of validity. Last, but not least, in Fig. 7(d) we 
plot the exact phase diagram for a system with Ny = 100, 
with the small boxed region corresponding to the region 
of validity for the multi-band approximation, indicating 
that only a very small part of the full phase diagram can 
be recovered using this approximation. 

We want to emphasis that although it is always pos¬ 
sible to construct an effective multi-band theory, except 
for the cases where only 1 or 2 channels are sufficient to 
describe the systems, which are indeed more tractable 
analytically, the exact model has in general quite an ad¬ 
vantage of use over the multi-band approximation. This 
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FIG. 7. (Color online) A comparison of the exact and ef¬ 
fective phase diagrams as a function of B and /a. We take 
A = 0.25^n* and a = In panels (a,b) we present 

a comparison of the full and multi-band effective topological 
phase diagram for Ny = 50 and n* = 6,8 respectively. The 
white (topologically trivial) and dark purple (topologically 
non-trivial) phases denote the regions where both calculations 
agree. The pink regions are predicted to be topological by the 
exact calculation and non-topological by the multi-band ap¬ 
proximation, while the light purple regions are predicted to 
be non-topological by the exact calculation and topological 
by the multi-band approximation. Panel (c) shows a compar¬ 
ison between the phase diagrams obtained by the multi-band 
approximation with n* = 10, for systems with Ny = 100, 
and Ny = 101 respectively. The small pink regions are the 
extra-topological regions corresponding to Ny = 101. Panel 
(d) shows the exact topological phase diagram for Ny = 100 
with the region where the multi-band theory is working (cor¬ 
responding to the results plotted in (c)) boxed by the small 
rectangle . 


is because it not only captures accurately the physics for 
a wider parameter range, but it also allows in general 
for a more straightforward calculation of the topological 
invariant. This is simply because the real space provides 
a simpler description of the transverse direction in the 
presence of both the spin-orbit coupling and of the open 
boundaries. 


VII. CONCLUSIONS 

We have studied the formation of Majorana bound 
states in square lattices with longitudinal and transver¬ 
sal Rashba, a Zeeman magnetic field, and in the presence 
of superconducting proximity. We have considered fully 
open systems, which we believe is the most relevant sit¬ 
uation for experiments. We have first calculated exactly 
analytically the topological phase diagrams for long quasi 
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ID wires using a topological invariant. We have shown 
that small regions of this phase diagram can be repro¬ 
duced by a multi-band approximation under certain con¬ 
ditions and for certain system parameters. We have then 
explored numerically the phase diagrams for fully open 
systems with arbitrary dimensions and geometries using 
the Majorana polarization. We have found that for very 
elongated systems the MP criterion recovers the phase di¬ 
agrams obtained analytically. For systems in which the 
length and the width of the system become comparable 
we find that more peculiar features arise: thus the low en¬ 
ergy states are no longer Majoranas but quasi-Majorana 
states, i.e. states that have a small but non-zero energy, 
as well as locally a perfect electron-hole overlap, but no 
global Majorana symmetry. For square systems no true 
Majorana states can form but the quasi-Majorana states 
form in the entire region corresponding to the bulk 2D 
topological phase. We have also examined the asymme¬ 
try between transversal and longitudinal Rashba, as well 
the transition between the periodic and open boundary 
conditions. 

It would be interesting to explore these quasi-Majorana 
states more thoroughly, in particular, given the likelihood 
to generate such low-energy states in experimental sys¬ 
tems, in would be useful to investigate their possible in¬ 
terest for braiding operations and quantum computation 
as well as their statistics. Also, it would be interesting 
to understand the relation between these states and the 
chiral Majorana states. 

It would be also interesting to generalize our work to 
quasi ID wires with a three-dimensional structure and 


with and more complicated lattice structures and tight- 
binding models,as well as to explore the modification of 
the direction of the magnetic field, to see if the discrete 
phase transitions in the phase diagram can reproduce the 
experimental data for In As and InSb wires. 
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Appendix A: Analytical expressions for the 
topological invariant for: A, = 2,3 

We give here the full expression for the topological in¬ 
variant for Ny = 2,3. Eor larger numbers of wires it 
becomes quickly unwieldy, though it can always be cal¬ 
culated using the method we have shown, even up to very 
large numbers of wires. Erom Eq. (9) in the main text 
we have the topoloigcal invariant 

(5 = sgndetH(fi)detH(f2). (Al) 

Eor Ny = 2 one finds, reinserting tx = ty = t for clarity 
and defining fk = —2tcosk — /i. 


^^(^,2) = 


/ —t /o,7r + B 

/o, 7 r + B 


— t 


0 

iay + A 


ICXqj 


0 


0 


— Idy + 

V—io^^ + A 0 


^ t —fo,7T + B 

“/o,7r B t 


This gives 


detn{T,,,) =B^- -B- + 2Ui^ + + A^ - + A^) 


for the expression for the topological invariant. 
Eor Ny = 3 one finds, again reinserting t. 




and hence 


0 

— t 

/o, 7 r + B 

0 


\ay 

A 

-t 

/o, 7 r B —t 

iay 
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—iay 

fo,7T + B 

—t 

0 

A 


my 

0 

0 

—iay 

A 

0 


t 

~/o,7r + 

lay 

A 

—iay 

t 

— f0,7T + B 

t 

A 

iay 

0 

-f0,7T + 

B 

t 

0 

- - 4ay + 3A^) 

> “ (/o,7r + 

(/ob 

+ 2/ot 

.(A^- 

2t^ - 2o? 


+ 8^2 

+ 3A" + /o% 

,(6A2 - 


5 



(A2) 


(A3) 


(A4) 


(A5) 


for the expression for the topological invariant. Appendix B: No transverse SO: = 0 

Clearly for weak transverse a the phase diagrams can 
be understood as a perturbation from the case of wires 
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which are transversally coupled only by hopping, without 
spin orbit interaction along the transversal directiond^ 
Here we provide some details corresponding to arrays of 
wires with = 0 , using the same notations as in the 
main text. An array of wires coupled only by a hopping 
term has been studied before. Here we simply rederive 
the topological invariant using the same method as in 
Sec. HI. After a Fourier transform with periodic bound¬ 
ary conditions (PBCs) imposed along x and OBCs along 
the system can be written as with 


(f{k) + B 


U{k) = 


where 


p * 

‘-'k 

-A 

0 


^k 

f{k) - B 0 
0 B - f{k) 

-A 


0 

-A 

n * 

^-k 


-f{k) - BJ 
(Bl) 


f{k) = -t{cos[kx] + cos[/c^]) - li , 


(B2) 


and 



= -i(asin[A:a,] (B3) 

is the spin-orbit term which vanishes at the two time- 
reversal-invariant (TRI) points F^i 2 ) = 0 , 7 r. Here k^ = 
27^1x1 Nx are the bulk momenta with n = 1 , 2 ,... and 
ky = 7rm/{Ny + 1) with m = l,2,...A^y are quantum 
numbers labeling a transverse degrees of freedom. 

The Hamiltonian for any ky at the TRI momenta can 
be written in block diagonal form as 


-Hk^Xi) 




_o ^ 

-nkATi)) 


where in this case 


nk^Ti) 


fB + f{ri,ky) -A \ 
I, -A B-m,ky)) 


(B4) 


(B5) 


The topological invariant for each ky is then, provided 

(a 7^ 0, 

Sk^ =sgn [defH{Ti,ky)defH{T2,ky)] 

= sgn { [B^ - A^ - f{0, ky)] [B^ - A2 - f^n, ky)]} . 

(B 6 ) 

When ^ = — 1 the system is topologically non-trivial and 
hosts Majorana states on the edges, for = 1 it is topo¬ 
logically trivial. 

Fig. 8 shows the topological phase diagram and the 
Majorana number as a function of B and y for several 
systems with ay = 0. Some regions of the phase dia¬ 
gram support a large number of MBS protected by chi¬ 
ral symmetry.In Fig. 9 we plot the number of edge 
states. Although the general pattern follows the form de¬ 
rived from the topological invariant, the quasi-ID wires 
show also a variety of unprotected edge states. 


FIG. 8 . (Color online) Topological phase diagram as a func¬ 
tion of B and y for (a,c) Ny = 3 and (b,d) Ny = 15, with 
transverse = 0 and A = 0.4. For (a) and (b) light red 
is topologically non-trivial and white is topologically trivial, 
where we have plotted the parity of the Z invariant which is 
a Z 2 invariant. Panels (c,d) show the number of MBS pairs, 
M. 
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FIG. 9. (Color online) Number of edge state pairs, A, as a 
function of B and y for a system with Ny =4, A = 0.4 and 
(a) ay = 0, (b) ay = 0.3. For ay = 0 many of these states 
are the protected MBS, which become gapped for 7 ^ 0. 


Appendix C: Asymmetric spin-orbit coupling: weak 
and strong coupling limits 

It is also straightforward to calculate the topological 
phase diagrams for the situation where the hopping and 
spin-orbit coupling is not of the same strength in the 
transverse and longitudinal directions. Consider Eq. ( 1 ) 
where we now allow ty 7 ^ tx and ay ^ ax- We will keep 
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FIG. 10. (Color online) Topological phase diagram as a 
function of B and /x for systems with tx = t, ty = 0.2t, 
o^x,y = 0.5tx,?/ and A = 0.4t. Light red is topologically non¬ 
trivial and white is topologically trivial. Plotted are phase 
diagrams for (a) 15 wires, and (b) 16 wires. 


the relative strengths of the spin orbit coupling in the 
two directions the same, ayja^ = tyltx, giving us just 
one more parameter to change: ty. We can calculate 
the topological invariant using the same methods as in 
Sec. Ill and will focus on two limits, weak transverse 
coupling ty ^ G, and strong transverse coupling ty ^ G. 

First let us condor the weak coupling limit ty <^tx. For 
an odd value of Ny it is clearly seen from Fig. 10(a) that 
the ID phase diagram is approached as ^ 0. In this 
limit the wires become almost decoupled in real space, 
each with a very similar expression for the bulk invariant. 
Thus in the region where all wires are topologically non¬ 
trivial there is an odd number of MBS at one end of the 
wire. When these states become coupled and hybridize 
one Major ana state must always survive. For an even 
number of wires this of course results in a topologically 
trivial region see Fig. 11(b), as the even number of MBS 
can all combine and destroy each other. 

Fig. ll(a,b) show the opposite limit where ty » tx. 
In this case the wires are all strongly coupled and the 
situation is similar to the isotropic case ty = tx. For a 
large ty the principle effect is to introduce a large energy 
scale such that the structure in the phase diagram is of 
the order of rather than t. Superficially the results 
are similar to Fig. 2. 



/i/t ^l|i 


FIG. 11. (Golor online) Topological phase diagram as a 
function of B and /x for systems with tx — t^ ty — 20t, 
OLx,y = t).btx,y and A = 0.4t. Light red is topologically non¬ 
trivial and white is topologically trivial. Plotted are phase 
diagrams for (a) 15 wires, and (b) 16 wires. 



012345012345012 

M M ju/t 


FIG. 12. (Golor online) Topological phase diagram as a func¬ 
tion of B and fi for a system with 3 wires with ax — Oiy — 
and A = 0.4. Light red is topologically non-trivial and white 
is topologically trivial. Plotted are phase diagrams for (a) 
t' — 0.3, (b) t' — 0.7, and (c) t' — 1 (the fully periodic nan¬ 
otube). Note these are all ID topological phase diagrams. 


for two examples. For a nanotube with an even Ny and 
t' = 1 the phase diagram is clearly superficially equiva¬ 
lent to the 2D phase diagram. Nonetheless the nanotube 
and square lattice are defined by different topological in¬ 
variants. As a function of t' the phase diagram is modi¬ 
fied by varying and merging gap closing lines at the TRI 
momenta. We can note that lines belonging to the two 
different TRI momenta can never merge as the longitu¬ 
dinal momentum is always a good quantum number. For 
t' = 1 as ^ oo the system resembles increasingly 
accurately a 2D system with 2 boundaries. The forma- 


Appendix D: From wires to cylinders 

The calculation in Sec. Ill remains valid if we consider 
attaching the outer wires 1 and Ny to form a periodic 
loop in the y direction creating a cylindrical geometry. 
We add to the Hamiltonian Eq. (1) a term 

Nx 

^' = E (-1 + + H.c., 

X=1 

(Dl) 

We can then calculate the ID topological invariant for 
ladders, cylinders, and for anything in between by vary¬ 
ing the strength of this connection, see Figs. 12 and 13 



M ju/t ju/t 

FIG. 13. (Golor online) Topological phase diagram as a func¬ 
tion of B and y for a system with 6 wires with ax = ay = 0.5 
and A = 0.4. Light red is topologically non-trivial and white 
is topologically trivial. The coupling t' has the same values 
as in Fig. 12. 
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FIG. 14. (Color online) MP for a system with 6 wires with 
OLx — OLy — 0.5 and A = 0.4. We show the transition from 
the MBS in an open system to a non-MBS in the periodic 
system. The state is destroyed by leaking along the weakly 
coupled edges. The coupling is t' — 0,0Aty,ty for (a,b,c) 
respectively. The system size is Nx = 60 and Ny = 6, see 
Fig. 13. 

tion of the 2D chiral Majoranas from the quasi ID edge 
states, as Ny is increased, has been thoroughly explored 
in Ref. 50. In the case under consideration here, as Ny 
is always finite, we remain always in the quasi-ID sce¬ 
nario and the edge states are modified by the coupling t' 
without forming a band. 

For an odd number of wires the nanotube possess only 


half the region of bulk non-trivial topology. This can be 
understood because for a periodic nanotube the topolog¬ 
ical invariant can be written in term of just four points 
corresponding to the TRI momenta for a 2D periodic 
square lattice. However for odd Ny only two of these 
exist in the actual system. 

By varying f we can switch from a quasi-ID ribbon to 
a quasi-ID cylinder and track the behavior of the lowest 
energy state using the MP. We have three interesting 
choices. Firstly if there is no topological phase transition 
and we start with a MBS then this state simply becomes 
translationally invariant along the transverse direction. 
No other change to the MP occurs. If, however we start 
with a MBS at t' = 0 and pass through a phase transition 
on the way to t' = ty then the MBS is destroyed, see 
Fig. 14. The lowest energy state gains a rotation in the 
MP along the transverse edge. Not that there are band- 
crossings as a function of f and hence the lowest energy 
state in Fig. 14(c) is not continuously modifiable into 
Fig. 14(a) by tuning t'. The MBS in Fig. 14(a) ends up 
as a bulk state for V = ty and the lowest energy state 
at t' = ty originates in an edge state at t' = 0. Finally 
we can start with no MBS at t' = 0 and pass through a 
phase transition on the way to t' = ty in which case a 
MBS is formed, this is simply the reverse of the previous 
case. 
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